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ABSTRACT 

When a self-gravitating disc is subject to irradiation, its propensity to fragmentation will be 
affected. The strength of self-gravitating disc stresses is expected to dictate disc fragmen- 
tation: as the strength of these torques typically decrease with increasing sound speed, it is 
reasonable to assume, to first-order, that disc fragmentation is suppressed when compared to 
the non-irradiated case, although previous work has shown that the details are complicated 
by the source of the irradiation. We expand on previous analysis of the Jeans mass inside 
spiral structures in self-gravitating discs, incorporating the effects of stellar irradiation and 
background irradiation. If irradiation is present, fragmentation is suppressed for marginally 
unstable discs at low accretion rates (compared to the no-irradiation case), but these lower ac- 
cretion rates correspond to higher mass discs. Fragmentation can still occur for high accretion 
rates, but is consequently suppressed at lower disc surface densities, and the subsequent Jeans 
mass is boosted. These results further bolster the consensus that, without subsequent fragment 
disruption or mass loss, the gravitational instability is more likely to form brown dwarfs and 
low-mass stars than gas giant planets. 
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1 INTRODUCTION 

Discs around both supermassive black holes and young stars are 
thought to undergo phases in which they are self-gravitating. It 
is possible that if these discs are sufficiently unstable while they 
are self-gravitating, they may fragment into bound objec ts, pro- 
■ viding a mode of planet formati on in protostellar discs ( Kuiper 
ll95ll:ICameronlll978l; Boss! 19971) or star formation in AG N discs 
jLevin & BeloborodovbooalNavakshin & Sunvaevll2005l) . 

The conditions under which a self-gravitating disc is expected 
to fragment have been extensivel y studied. Th e most important cri- 
terion is the Toomre parameter dToomrell 19641) 



1.7, 



(1) 



where c s is the local sound speed, K e is the local epicyclic fre- 
quency (equal to the angular velocity Q in Keplerian discs) and 
E is the disc surface density. This is a linear stability criterion - 
the critical value s given above ap ply to non-axisymmetric pertur- 
bations (see e.g. Durisen et al. 20071) . 

Once the disc becomes gravitationally unstable, spiral waves 
are excited in the disc, producing stresses and providing heating 
through shocks. This stress can be described as a pseudo-viscosity 
in some cases dShakura & Sunvaevl 1 19731 : iBalb us & Papaloizou 
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1 19991 : lLodato & Ricdl2004l iForgan et alj|201 ll) , allowing a sim- 
ple description of angular momentum transport in the disc. The 
stress produced by the spiral structure can therefore be described 
by a Shakura-Sunyaev turbulent viscosity parameter a, and semi- 
analytic models of self-gr avitating discs can be constructed a t 
low computational cost (e.g. lRice & Ar mitage 2009; Clarke 20091) . 
This pseudo- viscous approximation has bee n shown to fail if the 
disc is too massiv e or geometrically thick (Lodato & Ricell2005l : 
IForgan et aflEoI J), but it remains useful for moderately massive, 
geometrically thin discs. 

A balance can be struck between the shock heating from the 
spiral arms and the local radiative cooling. Discs whi ch achieve 
this b alance are described as being marginally stable dPaczvnskil 
Il978l) . This allows a relationship between the local cooling time 
normalised to the local angular frequency (/3 C = t coo iQ,), and a: 



(2) 



97(7 - l)/3c 



For fragmentation to be successful, a density perturbation pro- 
duced by a spiral arm should be able to grow until it is gravitation- 
ally bound. Therefore, the fragment should be able to cool effi- 
ciently to reduce pressure support due to thermal energy, and con- 
tinue collapsing. As cooling becomes more efficient in a marginally 
unstable disc (i.e. as /3 C decreases), a must increase to redress 
the balance. A consensus has developed in recent years that there 
is a maximum value of a that the disc can sustain before quasi- 
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stead y self-gravitating torques saturate: a cr it ~ 0.06 dGammiel 
1200 ll : iRice et all l2005h. T his is supported by numerical experi- 
ments dCossins et alj20 09) that confirm an inverse relationship be- 
tween fic and surface density perturbations: 



>oc 



rWc 



(3) 



Once the disc reaches a state such that a cannot increase to 
maintain thermal equilibrium, AE/E is typically of order unity 
(i.e. density perturbations become non-linear) and fragmenta- 
tion can occur. Rapid cooling also helps to prevent destruc- 
tive fragment-fragment collisions, as the initial fragment spac- 
ing w ill be such that collisions oc cur within a few orbital pe- 
riods (Shi osman &Begelmanlll989l) . Rapid cooling can then en- 
sure that colliding fragments do not become unbound as a re- 
sult of the collision. The combination of these criteria assure 
fragmentation can only occur in the outer regions of pro tostellar 

than r = 30 - 40 AU (Rafikov 



discs , typically at radii greater 
20051: IMatzner & Levirfl2005l: 



Whitworth & Stamatellos 2006; 
Bolev et alj 120061: (Stamatellos & Whitwort H |2008l; iForgan et all 



20091 : IClarkdl2009HVorobvov & BasuhOloHForgan & Ricefeoill) 
In the case of AGN discs, this critical r adius is approximately 
r = P ipe dLevin & Beloborodovl 120031 : iNavakshin eTaT] 120071 : 
lLevinll2007l : I Alexander et al.ll2008l) . 

This description, however, is incomplete. As fragmentation is 
consigned to the outer regions of self-gravitating discs, the local 
temperature is unlikely to be determined purely by gravitational 
stresses, and is likely to be governed by irradiation, either from the 
central object or from an external bath (e.g. envelope irradiation 
in embedded protostellar discs). In this regime, there is no longer 
a uniquely determined relation between a and /3 C , as irradiation 
provides an extra term to the disc's energy budget. 

It is somewhat intuitive to assume that adding extra heating 
will push the disc away from marginal instability by increasing 
the sound spee d, weakening self-gravity and therefore suppressing 
fragmentation. IMatzner & Levinl d2005l) demonstrate this using an 
analytic prescription for protostellar disc formation and evolution 
from Bonnor-Ebert Spheres, showing that fragmentation is typi- 
cally i nhibited for orbital periods lower than 20,000 years. Rafikov 
( 2009) also show that fragmentation is inhibited at large orbital pe- 
riod s when the disc a ccretion rate is low. 

ICai et al.l |2008) studied the effect of envelope irradiation in 
grid based hydrodynamic simulations with radiative transfer, show- 
ing that envelope irradiation suppresses the higher m spiral modes 
of the grav itational instability (a result to some extent predicted by 
IBoss| [2002. although they disagree on the role of convection in re- 
sisti ng this suppresion). 

Stamatellos & Whitworthl d2008l) compare stellar and back- 
ground irradiation using smoothed particle hydrodynamics (SPH) 
si mulations. The background irradiation behaviour is similar to that 
of lCai et all d2008h - stellar irradiation allows the outer disc regions 
(r > 30 AU) to cool sufficiently rapidly to fragment according to 
the minimum cooling time criterion, but they no longer satisfy the 
Toomre Q criterion. 

More recently. Ikratter & Murrav-Clavl d201 ll) have noted that 
weakening self-gravity is not generally an impediment to disc frag- 
mentation i f the full effects of mass infall are considered. As was 
shown by IStamatellos & Whitworthl |2008), the cooling time cri- 
terion can be satisfied quite easily in irradiated discs (although 
maintaining a low Q is difficult), and that low mass discs in the 
irradiation-dominated regime may be made more susceptible to 
fragmentation given the correct infall rate. This is consistent with 



ICai et alj j2008)'s finding that mild irradiation increases gravita- 
tional torques (and therefore increases the equilibrium mass infall 
rate). 

These results should be compared with l ocal shearing sheet 
simulations of irradiated discs dRice et al.l201 ll) . Under irradiation, 
there is no longer a fixed fragmentation boundary using either /3 C 
or q. Comparing to the non-irradiated case, irradiated discs can 
fragment for lower values of a, but still requires rapid cooling. 

With the straightforward criteria for fragmentation in the non- 
irradiated regime becoming less clear when irradiation is incorpo- 
rated, it is instructive to consider other, more generalised fragmen- 
tation criteria. We have developed such a criterio n, based on the 
local Jeans mass inside a spiral wave perturbation dForgan & Ricd 
1201 lh . If a fragment is to collapse and become bound, it must have 
a mass greater than the local Jeans mass. By measuring how the 
Jeans mass evolves with time, we can ascertain whether regions 
of a self-gravitating disc are becoming more or less susceptible 
to fragmentation. From this criterion, we are able to estimate not 
only when a disc fragments, but what the initial mass of the frag- 
ment should be, an important initial condition to models such as 
the "tidal downsizing" hypothesis of terrestrial and giant planet 
format ion dBolev et ai.l2010l : lMichael et al.l201ll ; lNavakshinl2"oTll 
l2010al lbl). 

This is somewhat similar to other work which has focused on 
the relationship between s piral arm width and the local Hill radius 
(Rogers & Wadsley 2011). These criteria reflect different compet- 
ing influences: the Jeans criterion compares self-gravity to local 
pressure forces, whereas the Hill criterion compares self-gravity to 
local shear. The advantage of both these criteria is that they are eas- 
ily generalised to cases difficult to describe or explain by traditional 
minimum cooling time/maximum stress criteria. 

In this paper, we apply the Jeans criterion to ID self- 
gravitating disc models where the effects of irradiation are ac- 
counted for. We compare these to models where irradiation is not 
present, to assess whether irradiation promotes or inhibits fragmen- 
tation. In section [2] we outline the Jeans criterion and how the ID 
disc models are constructed. In section [3] we describe and discuss 
the results from the models, and in section [4] we summarise the 
work. 



2 METHOD 

2.1 Calculating the Jeans Mass in irradiated spiral arms 

To calculate the Jeans mass inside a spiral density perturba t ion, w e 
adopt the same procedure as described in Forg an & Ricel (1201 lh ■ 
We assume the Jeans mass is spherical: 



Mi 
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Under the thin disc approximation 

Ppert = E per t/2_H~ = E ^1 H — — ^ /2H. 

If we assume the disc is marginally stable, we can obtain 

4^2tt 3 Q 1/2 c 2 s H 



M, = 



3G 



(4) 



(5) 



(6) 



In 



Forgan & Ricel J20I 1), we used the empirical result of 



ICossins et al. I d2009l) to estimate the fractional amplitude AE/E: 
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(7) 



As we are now considering irradiated discs, it is more appropria te 
to use the empirical relationship determined bv lRice et al.ld201lh : 



< 



>= 4.47Va. 



(8) 



For the disc to be in thermal equilibrium, the radiative cooling, vis- 
cous heating and irradiation heating must balance. The dimension- 
less cooling time, /3 C is modified thus: 



/3c = 



(9) 



ffflu(T*- 3^)7(7 -1)' 
where Ti rr represents the temperature of the local irradiation field. 

2.2 The Fragmentation Criterion 

To produce a generalised fragmentation criterion, we consider the 
timescale on which the local Jeans mass changes. We define: 



rv = ^o. 

Mj 



(10) 



If this quantity is small and negative, then the local Jeans mass 
decreases rapidly, and fragmentation becomes favourable. Equally, 
if this quantity is small and positive, the local Jeans mass increases 
rapidly, and fragmentation becomes unlikely. 

The value of T j estimates the number of local orbital periods 
that will elapse before fragmentation becomes likely (if at all). In a 
steady-state disc in perfect thermodynamic equilibrium, Vj — > oo, 
and fragmentation will never occur. Discs that are not in equilib- 
rium will assume non-infinite values of Vj, and fragmentation is 
either likely or unlikely. 

We assume — 5 < Tj < for fragmentation, but this is yet to 
be established empirically. This is equivalent to requiring a region 
of the disc to fragment in less than five orbital periods from the 
moment of measuring Tj. Substituting for H in equation {6J gives: 



M.j 



4%/27T 3 



Q 



3G 0(1 + 4.47^7i)' 



We can calculate Mj using the chain rule 

Mj 



8Mj . OMj ■ OMj . dMj 
Cs + -t^M + — a + -——< 



dc a 90 da OLc 

We will assume Q — a = = 0, and hence 



Mj = Mj3— 

Cs 



From equation (22) of iForgan & Ricel d201 ll) . we calculate 

c. = 1/2 / 9o7( 7 - 1)0 _ _1_ 
Cs \ 4 t coo 

and hence derive: 



r, 



(11) 



(12) 



(13) 



(14) 



(15) 



As in IForgan & Rice] ( l201lh . we assume that the maximum value 
for a saturates at some value q a at = 0-1, which is sli ghtly higher 
than the canonical value of 0.06 dGammidl200ll ; iRice et alj [2005). 
We take this value as a conservative estimate given current uncer- 
tainties regarding convergence of 3D simulat i ons of fragmenta tion 
dMeru & Batel |20 111 : lLodato & Clarke] |20 111 : IRice et alj|2012h . If 



the local value of a < a S at, then the system is in thermal equi- 
librium and Tj — s> oo. Once the local a reaches a sat , it is fixed at 
this value, while f3 c is allowed to decrease, reducing the local Jeans 
mass at a rate given by Tj. Comparing Tj to the critical range of 
values given earlier then determines whether a fragment shall be 
formed. 

We should note that we do not consider infall in this analysis, 
although we do assume the discs reach a steady state where the 
disc accretion rate M is the appropriate value to process material 
to maintain a constant disc mass (i.e. the local rate of change of 
surface density E = 0). From the above, we could suppose that a 
disc experiencing infall such that E > may give Tj < 0, and 
as such irradiation could encourage fragmentation provided that Q 
does not also increase. We leave investigation of this possibility to 
future work. 

It is unclear from this analysis whether irradiation should drive 
or stop fragmentation - consequently it is instructive to test this 
using semi-analytic disc models. 



2.3 Simple Irradiated Discs with Local Angular Momentum 
Transport 

As in Forgan & Rice ( l201lh . we co nstruct simple disc models in 
the same manner as Levin (2007) and Clarke (2009) to evaluate the 
dependence of the Jeans mass on disc parameters. We fix the disc's 
accretion rate 



M = 



37rac s E 
O 



(16) 



to be constant for all radii, and impose an outer disc radius. Assum- 
ing marginal instability (Q = 2) and thermal equilibrium then fixes 
the disc surface density profile. If thermal equilibrium demands 
a. > "sat, then a is set to a sa t, and Tj can then become a fi- 
nite, negative quantity. If, or when, —5 < Tj < 0, the disc model 
is assumed to fragment, and the local Jeans mass is calculated. 
We fix the star mass at 1M®, and consider four scenarios: 

(i) No irradiation (T irr (r) = OK), 

(ii) Envelope irradiation, Ti rr (r) = 1QK at all radii, 

(iii) Envelope irradiation, Ti rr (r) = 307^ at all radii 

(iv) Irradiation from the central star. 

For the third c ase, we use the radiation field expression given by 
lHavashilfl98l : 



T irr {r) = 280^ 



V 1M { 



IAU 



-1/2 



(17) 



For a sat = 0.1 and 7 = 5/3, we can actually determine the 
maximum value of /3 C at which fragmentation occurs (i.e. where 
T j = —5), which we find to be 



Pcrit < 2.2. 



(18) 



Note the conventional critical cooling time for this value of 7 is 
/3 C rit = 3. We have however assumed a value of a sa t around 60% 
larger than the usual value, so our slightly smaller value of /3 cr it 
compares sensibly to the typically used cooling time criteria. Al- 
tering the critical value of Tj instead would also allow /3 cr it to be 
increased to a more standard value, but this would permit discs to 
fragment on timescales longer than might be considered physical. 
In any case, the resulting fragment masses are only very weakly 
sensitive to these considerations. 
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3 RESULTS 
3.1 Disc Profiles 

Figure Q] shows 2D contours of the disc-to-star mass ratio q in all 
four scenarios, given the accretion rate and the disc outer radius. 
The differences between the no n-irradiated case (top left in Figure 
[]} and the models discussed in Forg an & Ricej d201 ll) are due to a 
different selection of Q (in this paper, we assume Q = 2, and in 
the previous work we assumed Q = 1), but the same broad features 
remain: a high disc-to-star mass ratio is required to maintain a self- 
gravitating disc with a modest accretion rate and reasonable outer 
disc radius. For example, in the non-irradiated disc, an accretion 
rate of 10 MQyr -1 and outer radius of r out = 30 au demands a 
disc mass of 0.119 M Q . 

At high accretion rates, the lifetime of material in the disc be- 
comes comparable to the orbital timescale. For this reason, we trun- 
cate the upper limit of the contours by demanding that 



M disc 2tt 

— : — > 5"FT> 
M fi 



(19) 



which essentially requires the disc to exist for at least five orbital 
periods at the given radius. Any disc model which does not satisfy 
this is discarded and not considered in the subsequent analysis. This 
explains the empty regions in the upper right portion of each plot. 
Comparing the non-irradiated case to the other three cases, 



we can see that the high M section of the parameter space (i.e., 
above M ~ lO^Moyr -1 ) remains similar. The effect of ir- 
radiation becomes more apparent at lower accretion rates, forc- 
ing the equilibrium disc structure to be more massive for a given 
M — r locus. For example, in the case of background irradiation at 
10K, the 0.214M© contour at r — 150 au moves downward from 
M = lO- 7 M yr- 1 to ~ 2 x lO _8 M yr- 1 . At T irr = 30A", this 
contour intersects the x-axis at < 100 au. Indeed, discs with radius 
150 au and accretion rates of lO _7 M0yr _1 have masses around 
0.357 M Q . 

In the stellar irradiation case, the change is less striking - a 
0.214Mq disc with radius 150 au is still formed at an accretion 
rate of ~ 2 x 10 _8 MQyr _1 , but the lower M region allows higher 
mass discs to have larger outer radii than the background irradiated 
cases. 



Adding extra mass to the outer regions of self-gravitating discs 
will in general enc o urage fragmentation dVorobvov & Basull20icl: 
Kratter et alj l2010l : iKratter & Murrav-Clavl 1201 ll; iForgan & Rice! 
2012b . so we might expect to see the Jeans criterion satisfied at 
lower accretion rates in the presence of irradiation. Conversely, in- 
creasing the disc sound speed will increase the typical fragment 
masses (if fragmentation can still be achieved). 
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Figure 2. 2D contour plots of the Jeans mass inside a spiral perturbation, as a function of the steady state 'pseudo- viscous' accretion rate M and radius in the 
disc r, for the four test cases considered in this work: without irradiation (top left), with background irradiation at 10K (top right), with background irradiation 
at 30K (bottom left), and with stellar irradiation (bottom right). Again note that the upper regions of the parameter space are excluded because the disc lifetimes 
are too short. Regions beneath the lowest contour are marginally unstable, but in thermal equilibrium, and do not fragment. 



3.2 Jeans masses 

Figure [2] shows the regions of parameter space in which the four 
scenarios can fragment, as well as the local Jeans mass at the point 
of fragmentation. 

The non-irradiated case is qualitatively similar to that found 
bv lForgan & Ricel feOllh . with the minimum Jeans mass also tak- 
ing a similar value of 2.66 Mj up . Self-gravitating discs will not 
fragment at low disc radii (except at very high, probably unphysi- 
cal accretion rates). There is a minimum accretion rate required for 
fragmentation even at large disc radii, corresponding in this case to 
around 7 x 10 M©yr _1 . By comparison with FigureQ] we can 
see this sets a minimum mass at which discs can fragment at around 
0.2A1q for a disc outer radius of 50 au. 

At high accretion rates, the regions of parameter space suscep- 
tible to fragmentation remain roughly similar across all four cases, 
but the irradiated discs show suppression of fragmentation at lower 
accretion rates. The 10K irradiation case shows fragmentation at 
accretion rates around a factor of 1.25 higher than the no irradiation 
case. This increases to a factor of 5 as the background temperature 
is increased to 30K. The stellar irradiation case displays similar be- 
haviour at low radii, but at around 70 au we see discs with low 
accretion rates possessing sufficiently low optical depths for the ir- 
radiation to heat the midplane so that Ti rr ~ T. While marginally 
unstable solutions still exist, fragmentation is suppressed at low 



accretion rates. At radii of ~ 120 au, accretion rates as high as 
10 _ Mqjt~ are required to produce a fragmenting disc. 

The minimum Jeans mass is typically larger for the three ir- 
radiated cases: 4.1 Mj up in the 1QK irradiation case, 15 Mj up in 
the 30 AT case and 1 1.2 Mj up for the stellar irradiation case. Also, 
the rate at which the Jeans mass increases with disc radius is higher 
than in the non-irradiated case. This further supports the creation 
of low m ass binary sta r systems as opposed to single-star planetary 
systems faarkell2009h . 

In short, irradiated discs tend to produce more massive frag- 
ments. There are two reasons for this: 



• The effective a generated by the instability is reduced in the 
presence of irradiation, and hence at a given M — r locus, the discs 
are more massive, and 

• The steady state temperatures in irradiated discs are higher, 
providing extra pressure support against gravitational collapse. 



Further , if we consider the expected isolation masses of these frag- 
ments jLissaued[l987h . then we expect irradiation to push frag- 
ments well into th e brown dwarf/low-mass star regime (see also 
lKratteretalj20ld) . 
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3.3 Surface Densities at Fragmentation 

We have established irradiation suppresses fragmentation at lower 
accretion rates, and boosts disc masses. If this is the case, the 
surface density at the point of fragmentation should typically be 
higher. In Figure|3] we again plot contours showing the leans mass 
of the fragments, but now study the E — r parameter space. As 
expected, lower mass fragments form at lower surface densities, 
and the fragment mass tends to increase with E. The disc lifetime 
criterion described in the previous section excludes part of the pa- 
rameter space, increasing the gap seen in the contours from approx- 
imately 50 au onwards. 

We see that irradiation does not affect the general trend at high 
E, but the low E fragmentation is suppressed. This is seen weakly 
in the background irradiation case, and much more strongly in the 
stellar irradiation case. Low surface density regions will generally 
be optically thin, and hence have their temperatures set by irradia- 
tion, remaining too warm to produce fragments. 



3.4 Cumulative Distributions 

Although a uniform parameter space in M and r is not what nature 
presents to us, it is instructive to see the resulting probability dis- 
tributions for fragment mass (and disc-to-star mass ratio) that can 
be inferred if we assume that each M — r combination in the pa- 
rameter space is equally likely. These can then be convolved with 
the true (M, r) distributions (once known) to predict the observed 
distributions of fragment properties. 

The left panel of Figure [4] shows the cumulative distribution 
of fragment masses found in this uniform parameter study. As 
low E regions do not fragment when irradiated, and disc masses 
are boosted, consequently the initial fragment masses are typically 
higher for the background and stellar irradiation cases (where the 
effect becomes significant for the stellar case). The peak of the mass 
distribution is shifted upward by around a factor of 2 (from around 
10 Mj U p to around 20 Mj up ) compared to the case without irradi- 
ation. 

With the surface density at fragmentation higher in the pres- 
ence of irradiation, and the disc masses typically higher for a given 
M — r locus, one might assume the typical disc-to-star mass ratio 
would itself be higher. In the right panel of Figure|4] we plot the cu- 
mulative distribution of disc-to-star mass ratios, assuming that all 
fragmenting discs extend to r = 150 au so as to make the compar- 
isons consistent. 

The resulting distributions depend on what fraction of the 
M — r parameter space permits fragmentation. From Figure |2] we 
can see the non-irradiated case permits fragmentation in the largest 
area of parameter space, and at the lowest accretion rates. As such, 
the disc masses required for fragmentation are the lowest, begin- 
ning at around 0.2 Mq (for the selected disc outer radius of 150 
au). The 10K irradiation case fragments at slightly higher M, and 
boosts disc mass slightly for the same accretion rate, and as such 
shows a slightly higher cumulative fraction. The 30K case frag- 
ments at noticeably higher accretion rates, as does the stellar irra- 
diation case. However, stellar irradiation suppresses fragmentation 
at large radii and low accretion rates, preventing lower disc masses 
from fragmenting. This does not occur for the 30K case, and as 
such its distribution is skewed towards even higher masses than the 
stellar irradiation case. 



4 CONCLUSIONS 

We have continued our investigation of the Jeans mass in spiral 
arm perturbations as a means of predic ting disc fragmentation and 
fragment masses (For gan & Ric el201ll) . extending the formalism to 
include the effects of irradiation, either from a constant background 
or from stellar irradiation. We construct simple one dimensional 
disc models, where the disc maintains thermal equilibrium between 
the heating due to self-gravitating spiral structures, irradiation and 
the radiative cooling. From this, we can discover the region of M — 
r parameter space in which fragmentation is expected to occur. We 
can also determine the expected fragment mass as a function of 
disc parameters. We compare three scenarios: non-irradiated self- 
gravitating discs, self-gravitating discs in a uniform temperature 
bath, and self-gravitating discs irradiated by the central object. 

Adding irradiation alters the equilibrium disc mass for a 
given accretion rate and outer radius, requiring higher mass self- 
gravitating discs to exist for these parameters. We find that in gen- 
eral, adding irradiation will suppress fragmentation at low accretion 
rates compared to the non-irradiated case. This is a similar conclu- 
sion to that drawn bv lRafikovl (2009). This is not surprising, as the 
analysis used here is similar, the only difference being the precise 
application of the maximum stress condition for fragmentation. 

Irradiated discs typically must sustain a higher accretion rate 
in order to fragment, and the local surface density at fragmentation 
is typically higher, as low surface density regions will not be gravi- 
tationally unstable due to the radiative heating. The fragment mass 
is typically boosted, and, as a result, disc fragmentation is more 
likely to form brown dwarfs and/or low mass stars than gas giant 
planets (unless some action of the disc or central star can prevent 
the fragments growing to their isolation mass). 

As this work investigates azimuthally averaged disc profiles, 
we are not able to address the possibili ty of irradiated stochas- 
tic fragmentation (cf |Paardekooperll20 1 21) . which suggests regions 
which experience a small density perturbation can be irradiated in 
such a fashion that the local cooling time is reduced, and fragmen- 
tation could subsequently occur. This may be further assisted by 
the op acity regime in which the perturbation resides dCossins et al.l 
l20lj. 

This stochastic fragmentation condition is equivalent to a lo- 
cal density perturbation having a reduced Jeans mass for a sufficient 
duration to allow self-gravitating collapse to produce a bound ob- 
ject. The key issue here is: can the local value of Q also remain 
low? As we fix Q when carrying out these models, we cannot an- 
swer this question. Future work should investigate the behaviour 
of the Jeans mass in numerical simulations, both with and without 
irradiation, to firstly establish how the Jeans mass evolves during a 
"realistic" fragmentation event, and also to test whether stochastic 
fragmentation is a sensible outcome of Jeans mass evolution. 

To some ex te nt, we therefore agree with 
iRratter & Murray- Ciavl J20T 1): irradiated discs that possess 
high accretion rates can still be driven to fragmentation, i.e. more 
massive discs are required with sufficient infall to maintain them, 
as well as high surface densities and temperatures sufficiently low 
that g ravitational instability can still arise. We would also agree 
with Stamatellos & Whitworthl (2008) that the last of these three 
conditions is likely to be the greatest obstacle to the fragmentation 
of irradiated self-gravitating discs. 
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Figure 3. 2D contour plots of the Jeans mass inside a spiral perturbation for (left) a disc irradiated at a background temperature of 10 K, and (right) irradiation 
from the central star, as a function of the disc surface density £ and radius in the disc r, for the four test cases considered in this work: without irradiation 
(top left), with background irradiation at 10K (top right), with background irradiation at 30K (bottom left), and with stellar irradiation (bottom right). Regions 
without contours do not possess a fragmenting solution. 
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Figure 4. Left:The cumulative fraction of fragment mass in the parameter space investigated, assuming all values of M and r ou t are equally likely to occur. 
Right: The cumulative fraction of disc-to-star mass ratio for all models which produced fragments (assuming all discs extended to the maximum radius of 150 
au). 
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